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ABSTRACT 

We present a new short-period brown dwarf candidate around the star TYC 1240- 
00945-1. This candidate was discovered in the first year of the Multi-object APO 
Radial Velocity Exoplanets Large-area Survey (MARVELS), which is part of the third 
phase of the Sloan Digital Sky Survey (SDSS-III), and we designate the brown dwarf 
as MARVELS- lb. MARVELS uses the technique of dispersed fixed-delay interferom- 
etery to simultaneously obtain radial velocity measurements for 60 objects per field 
using a single, custom-built instrument that is fiber fed from the SDSS 2.5-m tele- 
scope. From our 20 radial velocity measurements spread over a ~ 370 d time baseline, 
we derive a Keplerian orbital fit with semi-amplitude K = 2.533 ± 0.025 km s" 1 , period 
P = 5.8953 ±0.0004 d, and eccentricity consistent with circular. Independent follow- 
up radial velocity data confirm the orbit. Adopting a mass of 1.37 ±O.11M for the 
slightly evolved F9 host star, we infer that the companion has a minimum mass of 
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28. 0± l.5Mj up , a semimajor axis 0.071 ±0.002 AU assuming an edge-on orbit, and is 
probably tidally synchronized. We find no evidence for coherent instrinsic variability 
of the host star at the period of the companion at levels greater than a few millimag- 
nitudes. The companion has an a priori transit probability of ~ 14%. Although we 
find no evidence for transits, we cannot definitively rule them out for companion radii 



1. INTRODUCTION 



One of the first results to emerge from high-precision radial velocity (RV) surveys seek- 
ing substellar companions was the existence of a brown dwarf (BD) desert: a paucity of close 
(a < 5 AU) brown dwarf (l3M Jup < M < S0M Jup ) compani ons to solar-type stars, relative to more 
common stellar mass companions (IMarcy & Butlerll2000|) . Indeed, since they induce reflex ra- 
dial velocity semiamplitudes of many hundreds of meters per second, such brown dwarf com- 
panions have been w ithin the detection capabilities of these sur veys for over two deca des (e.g., 
Campbell et al.lll988|) . yet to date only a few dozen are known (|Reid & Metchevll2008|) . On the 
other hand, as instrumentation has subsequently improved, first Jovian, an d now terrestrial plan 



etary companions in similar orbits have been found in relative abundance (|Cumming et al.l 12008 



Mayor & Udryll2008l : [Mayor et al.1 12009). The brown dwarf mass regime represents an apparent 
minimum in the mass distribution of close companions to solar-type stars. 

Planetary companions are believed to form in circumstellar protoplanetary disks, whereas 
stellar companions are believed to form by concurrent collapse or fragmentation, so the brown 
dwarf desert is commonly interpreted as the gap between the largest mass objects that can be 
formed in disks, and the smallest mass clump that can collapse and/or fragment in the vicinity of a 
protostar. Such a gap was by no means guaranteed to exist, and is perhaps surprising. For example, 
numerous isolated BDs in star-forming regions have been found to possess proto planetary disks, 
akin to the disks of young stars ; suggesting that BDs form muc h as stars do (e.g., ICaballero et al 



2007 : Luhman & Muench)2008 : Scholz & Jayawardhanall2008 ). More generally, the mass function 
of isolated substellar obje cts in the field and clusters appears to be roughly flat in logM for masses 
down to at least ~20Mj up ( Luhman et alJboOO : Chabrierf2002), whereas i t is not clear what sets the 
upper limit for objects formed in protoplanetary disks ( e.g., IBossII2001| ; llda & Lirj|2004c iRafikov 
20051 : bodson-Robinson et al.lliooil : iKratter et aliboich . 



As such, details of the demographics of companions in the brown dwarf desert, including 
the aridity of the desert, the shape of the high-mass tail of the planetary companion mass func- 
tion and the low-mass tail of the stellar companion mass function, as well as how these prop- 
erties change with semimajor axis and primary mass, encode a wealth of information about the 



-4- 



poorly understood physics of star and planet formation. Additional processes such as tidal evo- 
lution and disk-planet m igration can also affect these properties (e.g., lArmitage & Bonnellll2002l : 
Matzner & Levinll2005|) . thus can be investigated via brown dwarf desert statistics. 



Unfortunately, despite its potential diagnostic power, and after more than twenty years of 
precision radial velocity surveys, very little is known about the brown dwarf desert, precisely be- 
cause brown dwarf companions are rare and so few such companions are known. The California 
& Carnegie Planet Search finds an occurrence rate of 0.7% ± 0.2% from their sample of ~ 1000 
target stars (Vogt et al. 2002, Patel et al. 2007), and the McDonald O bservatory Planet Search 
agrees, with a rat e of 0.8% ± 0.6% from a search sample of 250 stars ( Wittenmyer et al.l l2009). 
Gizis et al.l (1200 II) suggest th at brown d warfs might not be as rare at w ide separations (see also 
Metchev & Hillenbrand! |2004|) . although McCarthy & Zuckermanl (|2004|) find a low rate of occur- 
rence that is similar to that found for close separations. By extrapolating the mass functions of 
planets (on the low mass side) and stellar comp anions (on the high mass side) into the brown 
dwarf mass regime, Grether & Lineweavei ( 20061) find a mass of minimum occurrence (the driest 
part of the brown dwarf desert) at 3V^Mj up . They further suggest that the location of this mini- 
mum may scale with host star mass. For instance, the only known BD eclipsing binary is a "desert 
dweller", consisting of a ~ 60M hip BD with a ~35M Jup BD companion at a separation of 0.04 AU 



(IStassun et alJl2006L 120071) 



To make further progress on understanding the properties of the brown dwarf desert, a much 
larger sample of brown dwarf companions is needed. Furthermore, this larger sample must be 
drawn from a relatively uniform survey with a well-defined and homogeneous sample of primary 
target stars, so that the demographic properties of these companions can be reliably inferred. Given 
the occurrence rate of ~ 1 %, a survey of ~ 10000 stars is needed to detect of order 100 brown dwarf 
companions. Such an extensive survey would require a prohibitive amount of observing time with 
traditional echelle-based precision RV instruments, which can only target one object at a time. 
Furthermore, in many cases the RV precisions that can be achieved with these instruments are 
far better than are needed to detect brown dwarf companions, implying that this is not the most 
efficient application of these instruments. 



The Multi-object APO Radial Velocity Exoplanets Large-area Survey (MARVELS; iGe et al 

2008|) is a radial velocity survey of ~ 11000 stars (~ 10000 dwarfs and subgiants, plus ~ 1000 



giants) with 7.6 < V < 12 over time baselines of ~ 1.5 years, with a stated goal of < 30m s 1 
precision for the faintest stars. It operates as one of the bright-time survey compo nents of the Sloa n 
Digital Sky Survey (SDSS) III, following on the legacy of the original SDSS dYork et alJl2000l) . 
MARVELS uses the innovative instru men tal technique of a dispersed fixed-delay inter ferometer 
(DFDI; see, e.g.. brskine & Gel I2OO0I : lGel I2OO2I : be et all I2OO2T : Ivan Evken et all hold) in order 
to simultaneously observe 60 objects at a time over a three degree field of view with a single 
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instrument that is fiber fed from the SDSS 2.5-m Telescope (|Gunn et al.l I2006Q at Apache Point 
Observatory. The fibers are fed through an interferometer, and both interferometer output beams 
are sent through a spectrograph with a resolving power 12000, producing fringing spectra 
over the wavelength range ~500-570nm. Radial velocity information is imprinted in the phases 
of the fringes perpendicular to the dispersion axis of the spectrum due to a fixed variation in the 
interferometer delay along this direction. 

By virtue of the large number of target stars, as well as uniform selection criteria described 
below, MARVELS is well suited to probe for rare companions. MARVELS commenced operations 
with SDSS-III in Sep. 2008, and as of the end of the first year's data collection in Aug. 2009, had 
observed 780 stars with RV time series of more than 15 points. In this paper, we report the first 
MARVELS brown dwarf candidate, which we desig nate MARVELS - lb, detected in orbit around 
the starTYC 1240-00945-1 (Tycho-2 star catalogue: llfog et alJ EoOoh. 



2. OVERVIEW OF SDSS-III MARVELS TARGET SELECTION 

The overall scope of MARVELS will be described in detail in future papers; we present a 
brief outline here in order to provide the context for the field and target selection of the brown 
dwarf candidate. MARVELS has been designed with an RV precision goal of < 30 m s" 1 in order 
to be able to discover a sample of ~ 150 new exoplanets, within a homogeneous parent sample 
of searched stars. By choosing a sample of target stars using a limited number of well-defined 
selection criteria, our sample suffers from minimal and well-understood biases, and can increase 
the size of the largest statistically homogeneous exoplanet sample by a factor of a few over that 
currently available. 

MARVELS will run for six years during SDSS-III bright time, in a series of three cycles 
of self-contained two-year surveys. Each cycle will have a similar stellar target selection strategy, 
designed to give good survey coverage of FGK dwarfs and similar parent samples in each two-year 
cycle, although with different target fields. While in general this means only companions with up 
to ~ 1 .5 year periods will be detected, the advantage of this strategy is that we need not wait the full 
six years to gather enough epochs per star to detect companions. Also, this approach provides the 
opportunity to do major instrument upgrades at the end of each two-year cycle, without destroying 
the continuity of our RV measurements. 

In order to collect enough photons to achieve < 30 m s" 1 statistical RV precision, the stars 
we monitor must in general be brighter than V ~ 12, although the precision at a given magnitude 
depends somewhat on stellar parameters as well. For the 60 object multiplexing capability of the 
instrument during the first two years (Sep. 2008-Sep. 2010), we found most fields on the sky 
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were sufficiently rich to fill all the object fibers, so half of our fields were selected to include a 
reference star of 8 < V < 12 with a known RV signal (stable or planet-hosting). By recovering the 
RV of the reference star, we can verify that the instrument is sufficiently stable to detect planetary 
companions. To ensure survey observability across all right ascensions, the remaining fields were 
selected from ar eas with no referenc e stars. Finally, we also selected some fields in the Kepler 



survey footprint (Bo rucki et a l. 1997) in order to have the potential to leverage the exquisite Kepler 



photometry for any stars targeted by both surveys. 



In e ach individual field, we used the intersection of the GSC2.3 (|Spagna et al.l |2006) and 



2MASS dskrutskie et alJbooeh catalogs as our initial targets database, but because many of our 



fields are in the Galactic plane and contaminated by giants, we conducted a preselection program 
to identify and reject giants from the sample prior to beginning RV monitoring. First, we performed 
a rough cut in magnitude and color, accepting only stars with 9.0 < V < 13.0 and (J-K) > 0.29. 
The faint magnitude limit rejects stars too faint for the survey, and the bright magnitude limit keeps 
the dynamic range small enough to avoid saturating the preselection observations. The color cut 
eliminates most hot stars from consideration, since we cannot obtain sufficient RV precision to 
detect planetary companions on any star hotter than mid-F. Second, we took s pectral classification 



snapshots of the potential target stars using the SDSS double spectrographs (lUomoto et al. 1999) 



mounted on the SDSS 2.5-m Telescope, which have R ~ 1800 and cover the wavelength range 
390<A<910nm. 

The preselection observations were processed using t he SDSS two-dimensional and one- 



dimensional spectroscopic pipelines (|Stoughton et al.1 12002|) . The spectroscopic param eters r sff , 



logff, and [Fe/H] were derived using the SEGUE Stellar Parameter Pipeline (SSPP; Lee et al. 
2008). Each spectrum was manually inspected to validate the parameters and to identify obvi- 
ous binaries and emission line objects. 

The final 60 targets for each field were selected using the following method. First, we only 
consider stars with r e g < 6250 K. We dedicated 6 of the targets for observing giants, and identified 
the brightest available dwarfs and subgiants from 7.6 < V < 13.0 to fill the other 54 targets, where 
dwarfs and subgiants are defined as having logg > 3.0. For 7.6 < V < 9.0, we selected the targets 
for observation by conducting a literature search to reject known variable stars, and used a reduced 
proper motion (RPM) diagram to classify them as giants or dwarfs. While we prefer to pick 
bright dwarfs, in practice this bright magnitude range is dominated by giants, and therefore the 
MARVELS giant sample is typically drawn from the bright magnitude bin. For 9.0 < V < 13.0, we 
ranked the stars by V magnitude then picked the ~54 highest-ranked stars, although to avoid the 
survey being dominated by F-stars, we cap the number of stars with 5800 K< r eff < 6250 K at no 
more than 24 out of 60. In practice, this combination of criteria usually completes our 60 target 
selection without going fainter than V ~ 11.5- 12.0. We do not impose selections based on the 
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ages, activity levels, or metallicities of the stars. 

We have recently learned that the original version of the SSPP code that we used for our target 
selection tends to overestimate logg, particularly for cool temperatures of T eff <5000K. While we 
are working on improved methods to better discriminate between dwarfs and giants for targeting 
in future survey cycles, our target sample for the first two year survey cycle is likely to have more 
giants than we desired; we estimate that up to 30% of targets in this sample could be giants due to 
the bias in the SSPP results. Note that we do not use the primary properties of TYC 1240-00945-1 
derived from the SSPP in our subsequent analysis; we rely on the more accurate determinations 
from the detailed analysis of high-resolution spectra as described in §[6l We only describe the 
SSPP target selection method here because our MARVELS targets for years 1 and 2 (including 
TYC 1240-00945-1) have been selected based on the SSPP results. 



3. OBSERVATIONS AND PROCESSING 

3.1. Primary Survey Observations with SDSS 

TYC 1240-00945-1 was part of the first two-year cycle of the SDSS-IH MARVELS planet 
search program described above. This target was selected for radial velocity monitoring using the 
preselection methodology and instrumentation described in $2l In preselection observations for 
this star's field taken on Sep. 19, 2008, we obtained a series of five 7 s and five 12 s exposures 
of the target field, plus flat and arc lamp calibration exposures before and after this series. From 
preselection, the star appeared to be a late F-dwarf (but see further details in §(6} which suggest it 
is starting to evolve into a subgiant) suitable for inclusion in the MARVELS RV monitoring. 

Our discovery radial velocity observations were taken using the SDSS 2.5-m Telescope at 
Apache Point O bservatory coupled to the MARVELS instrument, a 60 object fiber-fed DFDI 



(|Ge et al.ll2009|) . Our two-output interferometer produces two fringing spectra ("beams") per ob- 
ject, over wavelengths ~500-570nm, with resolving power R ~ 12000. The instrument is envi- 
ronmentally stabilized such that no iodine cell is needed in the stellar beam path, and instrument 
drift calibrations are simply taken before and after each stellar exposure. TYC 1240-00945-1 was 
observed at 20 epochs from Nov. 7, 2008 to Nov. 11, 2009, as listed in Tabled] Exposures were 50 
min., yielding an average of 500 photons per CCD pixel on each 4k x 20 pixel fringing spectrum. 
The RV signal on TYC 1240-00945-1 was easily detected by eye in the RV curves from the first 
year of MARVELS. 

MARVELS RVs are differential measurements, based on the shift of the fringing spectrum 
relative to a template epoch. The RVs were derived from our 20 fringing spectrum observations 
using the preliminary version of our MARVELS DFDI pipeline, which is based on software from 



-8- 



earlier DFDI prototype instruments (e.g., Ge et al. 2006). We provide here a brief outline of the 
mechanics of the MARVELS-specific pipeline, but leave a full description to future techniques 
papers to be written on the overall performance of the MARVELS hardware and survey. 

After performing standard multi-object spectroscopic preprocessing on each frame such as 
bias subtraction, flatfielding, and trimming out individual spectra, we proceed to straighten slanted 
spectral lines, straighten tilted traces, and divide out uneven slit illumination to produce clean 
images ready for analysis. To remove a faint pattern of background fringes caused by the interfer- 
ometer, we apply a low-pass filter, which leaves the fringes on stellar lines visible. The pipeline 
seeks to measure the epoch-to-epoch shift in the two-dimensional fringing spectrum (i.e., a spec- 
trum with sinusoidal modulations along the slit direction). The shift induced by a stellar radial 
velocity change comprises two orthogonal components. The first component, a small shift of the 
stellar absorption lines along the wavelength axis, is the shift that conventional Doppler planet 
search instrumentation seeks to measure. The second component, a shift of the fringes on each ab- 
sorption line along the spectrograph slit axis, is linearly proportional to the shift in the wavelength 
axis, but is amplified to a factor of a few times larger, and therefore provides most of the statistical 
leverage in our velocity measurement. At any given wavelength, the fringe shift is related to the 
radial velocity by a multiplicative factor derived from measurements of the interferometer in the 
lab before commissioning the instrument. 

We use x 2 minimization to determine the best-fit velocity shift for each epoch, relative to a 
template spectrum chosen to be the brightest one from the epochs that were observed. Specifi- 
cally, we determine the best-fit velocity shift that minimizes the shift of the spectrum along the 
wavelength and slit axes, relative to the template spectrum. We also account for the barycentric 
correction during the RV extraction routine, ensuring that the x 2 minimizer does not need to search 
as far in velocity space as it would if the Earth's motion were not removed. Wavelength and slit 
axis shifts between exposures induced by the instrument drift were measured from fringing spec- 
tra of a stable calibration source (a tungsten lamp shining through a temperature- stabilized I 2 gas 
cell) taken before and after each stellar exposure, and the RV corrections due to these shifts are 
subtracted from each spectrum. Because the epoch for the instrument drift RV zero-point differs 
from the epoch for the stellar RV zero-point, none of these differential RVs will have a value of 
exactly zero. 

Because the interferometer splits the beam of each star, we record two separate spectra of 
each star on the CCD, and measure the RV from each of these spectra independently. We shall dif- 
ferentiate between these two simultaneously observed RV curves by using the labels "beaml" and 
"beam2." Although not all of the potential sources of systematic error would cause differences be- 
tween the two beams' RV curves, comparison of the two beams does provide a partial consistency 
check of the quality of the data and the reduction pipeline. 
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3.2. Photometric observations 



In order to check for intrinsic photometric variability indicating activity, as well as search for 
transits of the companion, we extracted the photometric time series data of TYC 1 240-00945-1 



obtained by the Kilodeg ree Extremely Little Telescope (KELT) North transit survey (Pe pper et al. 



20071 : ISiverd et al.ll2009l) . KELT consists of a 42 mm lens imaging a 26° x 26° field of view onto a 
4k x 4k CCD. KELT uses a red-pass filter with a 50% transmission point at 490 nm, which, when 
folded with the CCD response, yields an effective bandpass similar to R, but broader. 



The KELT data were processed as described in detail in Fle ming et al.1 (|2010 ). Brie fly, after 
flat-fi elding, relative photometry was extracted using the ISIS image subtrac tion package ( Alard & Lupton 
1998), combined with point-spread fitting photometry using DAOPHOT (|Stetsonl ll987). We re- 



duced flieJcvelof_systematics present in the light curve by applying the Trend Filtering Algorithm 
(TFA; lKovacs et alj2005|) . A few additional outlying measurements were removed before and after 



application of TFA. Raw uncertainties on the individual points were scaled t o force an ensemble o f 
stars near the target to have a modal x 2 /dof of unity for a constant fit. As in Fleming et al.l (|2010h . 
the target's x 2 /dof was still not unity after this adjustment based on the ensemble, so we further 
scaled the target's error bars by a small amount (~ 10%) to force x 2 /dof=l. The final KELT light 
curve has 5036 data points taken between Nov. 15, 2006 and Jan. 17, 2010, with typical relative 
photometric precision of ~ 1%. 



The Tycho catalog magnitudes (IHOg et al.l |2000() of our targets are generally unreliable at 
V > 1 1 . In particular, we find the error bars can sometimes be underestimated at the level of several 
tenths of a magnitude, resulting in colors that do not agree with spectroscopically-determined 
values of r eff . Therefore, we obtained absolute photometry to supersede and supplement the catalog 
colors. TYC 1240-00945-1 was observed in BV under photometric conditions by the privately- 
owned Hereford Arizona Ob servatory (HAO ) 11 -inch on Jul. 29 and 31, 2009, together with 
a program of standards from lLandoltl (|l992|) . It was observed again by this telescope in g'r'i' 
under photometric conditions o n Jan. 12, 2010, to gether with a program of Landolt standards that 
had u'g'r'i'z' calibrations from ISmith et al.1 (I2002|) . This telescope is equipped with a 1.5k x lk 
CCD with a plate scale of 0.81" per pixel. For each program night, the s tandard star instrumental 
magnitudes were fit with a generic photometric equation (see lGaryil2010i for more information on 
calibration procedures at HAO), and the resulting fit used to calculate the apparent magnitudes of 
TYC 1240-00945-1; typical standard star residuals relative to the fit were 0.01-0.02 magnitudes. 
The resulting calibrated BV g'r'i' are provided in Tabled Magnitudes in R r I r were estim ated from 
the measured g'r'i' by using the transformation equations tabulated in lSmith et al.l (|2002|) : the R C I C 
estimates are also provided in Tabled 
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3.3. Spectral Classification of Host Star 



In pursuit of a more detailed spectral classification of our candidate than is possible from our 
low-resolution SDSS spectrograph preselection observations, optical (~3600- 10000 A) spectra 
of TYC 1240-00945-1 were obtained on Nov. 2, 2009 with the Apache Point Observatory 3.5- 



m telescope and ARC Echelle Spectrograph (ARCES; iWang et al.l |2003|) . We used the default 
l."6 x 3"2 slit to obtain two moderate resolution (R ~ 31500) spectra with signal-to-noise ratio 
(S/N) of ~ 160 per 1-D extracted pixel at 6500 A. We extracted our APO classification spectra 
to 1-D using standard IRAF techniques and wavelength calibrated using ThAr lamp exposures 
obtained immediately after each science exposure. 

We also used the high resolution (/?=48000) spectrograph FEROS (|Kaufer et al. Il999h mounted 
at the MPG/ESO 2.2-m telescope in La Silla to obtain spectra of TYC 1240-00945-1. Two spectra, 
exposed for 3600 s and 4200 s respectively, were obtained in the wavelength interval 3500-9000 A, 
yielding a S/N ~340 per 1-D extracted pixel around 6600 A. These spectra were analyzed using 
the online FEROS Data Reduction System (DRS) and the standard calibration p lan, where bias. 



flat-fi eld and wavelength calibration lamp frames are observed in the afternoon. Idel Peloso et al 



(|2005|) checked the performance of the DRS by comparing the equivalent widths derived from solar 
spectra (observation s of reflected sunlight from Ganymede) with those from the Solar Flux Atlas 
(|Kurucz et al.lll984|) . They found that the two sets of measurements are strongly correlated, with 
a correlation coefficient of R = 0.994 and a standard deviation of 2.9 mA. The FEROS pipeline 
equivalent widths may thus be regarded as very robust. Furthermore, as the wavelength shift be- 
tween the two observed spectra was found to be negligible (13.1m s" 1 ), the two spectra were 
simply combined and shifted to the rest wavelength. 



3.4. Radial Velocity Follow-up 



To confirm this first substellar companion from MARVELS, as well as ascertain the qual- 
ity of the radial velocities obtained with the MARVELS instrument relative to those measured 
using convention al echelle spectrograph technology, we used the Hig h Resolution Spec trograph 
(HRS; lTulllll998l) mounted on the 9-m Hobby-Eberly Telescope (HET; iRamsev et aO 1998b to ob- 
tain additional precision RV measurements o f TYC 1240-00945-1. The candidate was observed in 
queue- scheduled mode (IShetrone et al. with a Director's Discretionary Time allocation espe- 
cially for this candidate, allowing for high-priority confirmation using just a few short (~ 15-min.) 
exposures spread over several nights. Nine measurements were taken in Dec. 2009 using an iodine 
cell for wavelength calibration, as well as one iodine-free template observation. All spectra were 
taken with the 316 lines mm" 1 grating with a central wavelength of 593.6 nm, leading to a resolv- 
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ing power R ~ 60000 and wavelength coverage 409 < A < 782 nm. Differential RVs were extracted 
from the HET spectra using a preliminary version of a new precise Doppl er reduction pipeline 



(kindly provided by Debra Fischer) based on the principles outlined in iButler et alj (|1996|) . This 
version of the pipeline was not yet optimized for the HRS fiber-fed spectrograph, and in particular 
used an instrumental profile description more appropriate for the slit-fed Hamilton spectrograph 
at Lick Observatory. As a result, systematic errors in the radial velocities presented here are high, 
and do not reflect the full capabilities of either the iodine technique or the HRS. The final measured 
radial velocities are given in Table [3] 

In addition, absolute radial velocities were obtained from the SMARTS 1.5-m telescope at 
CTIO. The target was observed 9 times from Aug.-Dec. 2009 using the echelle spectrograph with 
no iodine cell, yielding R ~ 42000 and wavelength coverage 402 < A < 730 nm. Each observation 
spanned 30 minutes of total exposure time, subdivided into three 10-min. exposures for cosmic 
ray removal. RVs were extracted using an IDL based pipeline written by F. Walter and adapted by 
K. Stassun. The individual exposures were bias-subtracted, flat-field corrected using quartz lamp 
flats, and wavelength calibrated using ThAr lamp exposures bracketing the science exposures. Typ- 
ically, 35 good echelle orders spanning 4800-7100 A, with a resolving power R ~ 42000, were ex- 
tracted from each observation, with a typical S/N ~ 30 per resolution element. Absolute RVs were 
measured via cross-correlati on against an early-K giant radial velocity standard star, HD 223807, 



selected from the catalog of iNidever et al.l (I2002|) . which was observed with the same instrument 



with S/N ~ 100. For each observation of TYC 1240-00945-1, cross correlation was performed 
order by order against the template, and the resulting 35 RV measurements from the individual 
orders were subjected to a sigma clipping based on the median absolute deviation. After clipping, 
we typically were left with RV measurements from 20-25 orders, which were averaged for the 
final RV measurement at that epoch. The measured absolute RVs are given in Table |4l We also 
applied this procedure to determine the RVs for six observations of the RV standard star obtained 
over the same time period; we found the root mean square (RMS) scatter in the standard star's RV 
measurements was 70 m s" 1 , which we take as the current precision limit of RVs obtained with the 
SMARTS 1.5-m echelle, without the Iodine cell and with the current preliminary pipeline. Note 
that, although the radial velocity standard is of a different spectral type than TYC 1240-00945-1, 
we expect that the systematic error that this mismatch produces will manifest itself primarily as an 
offset of ~ 1 - 2 km s" 1 added to all the absolute RV measurements, with a much lesser effect on 
the values of the RVs relative to each other. 
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4. RADIAL VELOCITY ANALYSIS AND KEPLERIAN ORBITAL SOLUTION 

4.1. MARVELS RADIAL VELOCITY DATA 

In Table [Q we present the 20 radial velocities measured by the MARVELS instrument, and 
we show the RV curve as a function of time in Fig. [T] Both beams are shown, and even though 
the error bars plotted in Fig. Q] are photon-only and do not account for systematics (our procedure 
to determine more realistic error bars follows below), several of the beamwise pairs nonetheless 
agree within their error bars. 

The MARVELS pipeline is still under development, and we find the RV scatter for other stars 
in the same field as TYC 1240-00945-1 (as well as for stars in other fields) is on average 2-3 times 
larger than the photon noise, on timescales greater than a month; presumably, most stars observed 
are not astrophysically variable at this level, indicating that the excess scatter is due to systematics 
(note we will discuss the expected RV jitter for TYC 1240-00945-1 in Section [6~4l since we need 
to determine the stellar properties first before searching for cases of similar stars in the literature). 
During pipeline development, we have examined the morphology of the RV residuals in the cases 
of several reference stars with known RV curves (either stable or planet-bearing) and found the 
systematic errors typically manifest in the form of month-to-month offsets at the level of tens of 
m s -1 , such that the RV data within any individual month fits the known RV curve much better than 
over multiple months. The offsets are often the same in direction and magnitude for both beams. 
These systematic errors may be due to imperfections in the detailed preprocessing of the images, 
because we do not see these systematics at the same level when analyzing simulated stellar data 
free of real- world image distortions. Since the exact factor by which the scatter exceeds the photon 
noise varies from star to star, we have decided to determine the excess scatter for the candidate at 
hand, to ensure that it falls in the typical range seen for other stars, and so is not responsible for 
the RV signal which we have interpreted as due to a companion. 

Our procedure for estimating the magnitude of the systematic errors in the RV curve is as 
follows. We assume that the systematic errors can be well-modelled by applying a simple constant 
multiplicative scaling to the uncertainties derived from the photon noise alone. We choose to use 
a multiplicative scaling of the error bars instead of adding a systematic error in quadrature to the 
statistical error bars because during pipeline development, we found the increase in RV scatter 
above the photon noise level is larger for fainter stars than brighter stars, so adding systematic 
error in quadrature would not be able to capture the overall form of the extra error as a function 
of signal-to-noise ratio. We designate this multiplicative scaling factor the "quality factor" Q. We 
estimate Q by performing a Keplerian fit to the RV dataset (with the raw pipeline photon-noise 
uncertainties), allowing for a linear trend with time. We then find the value of Q such that the 
X 2 /dof of the best-fit is equal to unity. 
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Following this error bar growth procedure, we found that the MARVELS RVs for TYC 1240- 
00945-1 were affected by systematics at levels of <2beami = 2.21 and <2beam2 = 3.63 for the two 
beams, respectively. Multiplying the statistical error bars by Q, we get a median scaled error bar 
of 92 m s" 1 for beam 1 and 151 m s" 1 for beam 2. After scaling the error bars, we performed a 
joint fit to beam 1 and beam 2 to provide a stronger constraint than a fit to either beam alone. The 
joint fit allows for different slopes and offsets between the two beams. This model is required 
because the two beams traveled through different parts of the instrument, and most importantly, 
experienced different optical path delays inside the interferometer (recall from Section 13.11 that 
there is a multiplicative factor that transforms fringe shift into radial velocity- this factor depends 
on the delay). The parameters of this final joint MARVELS orbital fit are given in Table [5] below, 
and the fit is overplotted with the data in Fig. [TJ The uncertaintie s were determined using the 



Markov Chain Monte Carlo (MCMC) method (see, e.g.. lFordll20 06). Note the time is referenced 
to the time of inferior conjunction (i.e., the expected time of transit if the system is nearly edge- 
on), and i s given as the Bary c entric Julian Date (BJD) in the Barycentric Dynamical Time (TDB) 



standard ((Eastman et al.ll2010|) . 



The Q values for the two beams are consistent with that of a typical constant star's Q, ~2-3. 
We also checked the brighter planet-bearing reference st ar HIP 14810, whic h was observed on the 
same plate at the same time. Using the known RV model ( Wright et al. l2009h . we find this reference 



star has Qbeami = 5.29 and <2beam2 = 4.36, with median statistical error bars of 9.1 m s" 1 for beam 
1 and 10.4 m s" 1 for beam 2. The higher Q for the brighter star is not an especially surprising 
result, since systematic noise sources that are independent of photon counts contribute a higher 
fraction of the total error when photon noise is small. Fig. Q] shows the residuals of HIP 14810 
relative to the model curve, on the same scale as the RV residuals of TYC 1240-00945-1. These 
residuals demonstrate that we can recover the RV curve of a known planet-bearing star to a level at 
least as good as our TYC 1240-00945-1 fit. Hence, the level of systematic uncertainty we find for 
TYC 1240-00945-1 is not unusual for its field, and that level is small compared to the amplitude 
of RV variability we find for TYC 1240-00945-1 and attribute to a companion- MARVELS- lb. 



4.2. HET AND SMARTS RADIAL VELOCITY DATA 

To further confirm that RV variability is indeed due to a companion, as well as to confirm 
the basic parameters of the orbital fit, we compared the RV observations obtained from HET and 
SMARTS to those obtained by the MARVELS instrument. We found these RV data do verify the 
variability and periodicity, but the follow-up data sets comprise insufficient high quality data points 
to provide much additional refinement of the orbital fit parameters on top of the discovery data. 

We first treated each RV dataset independently, computing a separate orbital fit and estimating 
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Q for the dataset using the procedure described above in 34.11 This gives the minimal error bars 
that would be consistent with any Keplerian orbital solution. We use these separate fits only for 
estimating the HET and SMARTS total error bars. 

For the HET data we find Qhet = 15.3, which is high, but expected due to the preliminary 
nature of the pipeline used to reduce the data (see 33.41) . We subtracted the RV model based on 
the MARVELS fit from the HET points and found that the residuals could be fit by a straight line 
(slope and offset) with % 2 = 7.9 and 7 degrees of freedom. Under the assumption that the errors 
are independent and normally distributed, this corresponds to a 33.4% probability of happening by 
chance, so there is no evidence to reject the hypothesis that the HET RVs are consistent with the 
MARVELS orbital fit. 

For the SMARTS data we find Qsmarts = 1-50. We subtracted the RV model based on the 
MARVELS fit from the SMARTS points and found that the residuals could be fit by a straight 
line (slope and offset) with x 2 = 20.4 and 6 degrees of freedom. Again assuming independent and 
normally distributed errors, this has a 0.23% probability of happening by chance, so there is strong 
evidence to reject the hypothesis that the SMARTS RVs are consistent with the MARVELS orbital 
fit. However, given that the HET and MARVELS RVs agree, we expect this discrepancy with the 
SMARTS data merely reflects evidence for unidentified systematics in the SMARTS data, which 
is not surprising, given the preliminary nature of the reduction of the SMARTS data (see 33.41) . 

The four RV data sets are shown in Fig. [2} phase-folded to the fitted period and phase (as 
determined from the fit to the MARVELS data alone). This visually demonstrates the conclusion 
that the HET and SMARTS RV data confirm both the amplitude and phase of the variability. We 
then tried an orbital fit to all three telescopes' data sets jointly, applying the same method that was 
used to jointly fit MARVELS beams 1 and 2, but now expanded to accommodate four RV data sets. 
We found that the new period and amplitude derived, using all the data sets combined, matched 
the values adopted in Table[5]to within the la uncertainties, and furthermore, that the uncertainties 
themselves matched to within ~10%. 



5. MONITORING FOR PHOTOMETRIC VARIABILITY 



The KELT data for TYC 1240-00945-1 are displayed in Fig. [31 and show no evidence for 
variability. The final weighted RMS i s 0.92%. A weighted Lomb-Scargle periodogram with float- 
ing mean (|Lomblll976c IScarglelll982|) yields no significant peaks for periods of 1 - 10 d, and in 
particular no evidence for any periodic variability near the period of the companion or the first 
harmonic. The improvement in x 2 for a sinusoidal fit at the period of the companion is only ~ 0. 1 
relative to a constant flux fit. 
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Fig.|4]shows the KELT light curve phased to the best-fit period of the companion (5.8953 d), 
as well as the phased light curve binned every 0.04 in phase (roughly the expected transit duration 
for a mid-latitude transit). The RMS of the binned curve is 0.059%, with a x 2 /dof of 0.85. This is 
consistent with no correlated (red) noise at the level of the RMS, since with an average of ~200 
data points per phase bin, one would expect a factor of ~ 15 improvement for the binned RMS 
compared to the unbinned RMS. We can also place an upper limit of 0.050% on the maximum 
light curve variability at a period half that of the period from the RV orbital fit (at A% 2 = 9), but 
this limit is insufficient to detect the exp ected amou nt of ellipsoidal variability for this candidate 



system. Using the equation in Table 2 of IPfahl et al.1 (|2008l) . we calculate the ellipsoidal variation 



would only be 0.0019% in amplitude. Note the methods we use to calculate the physical parameters 
for the star and companion used in the equations in this section will be explained later, in §16.31 



We possess an ephemeris from the RV orbital fit to search for companion transits at the ex- 
pected time. However, prior to our exposition of the Monte Carlo analysis using the RV infor- 
mation, let us first consider approximately what S/N to expect, calculated under the simplify- 
ing assumption of a random ephemeris (allowing us to write an analytic expression for the S/N). 
Based on the semimajor axis of a = 0.071 AU for an edge-on system, the a priori transit probabil- 
ity for the companion is fairly high, R*/a = 14.4%. The expected duration of a central transit is 
~ R*P I (na) = 6.49 hours, and the expected depth is 5~(r/_/?*) 2 = 0.2lS%(r / R Jup ) 2 , where r is the 
radius of the companion. Using these values, the expected S/N of a transit in the KELT data can 
be estimated, 

S/N~N^(^X /2 t^ 3 .5(-^-) 2 (1) 

where /V = 5036 is the number of data points, and cr~ 1% is the typical uncertainty. Thus the 
detection of a transit using KELT data is challenging if the radius of the companion r<R Jup , as 
is expected based o n the likely age of the star ( §17.11) and the minimum mass of the companion 



(Baraffe et al. 2003) 



in 



Detailed lim i ts on transits are produced by using the same Monte Carlo analysis as described 



Fleming et al.1 (|2010|) to incorporate our transit ephemeris from the RV data. Briefly, we use 
the distribution of companion periods and expected transit times from the MCMC chain derived 
from the fit to the MARVELS RV data (§|4J} to predict a distribution of transit times in the KELT 
data. For each link of the MCMC chain, we consider the uncertainty in the inferred radius of the 
primary due to the uncertainties in the spectroscopically measured r eff , logg, and [Fe/H] (see §16.31) . 
and we also consider a uniform range of transit impact parameters. For a given assumed radius for 
the companion, f or ea ch link we can then compute the expected transit curve using the routines of 



Mandel & Agoll (|2002|) . which are fit to the KELT dataset, computing the difference in \ 2 relative 
to a constant flux fit to the data. This is repeated for each link in the Markov chain, as well as for 
a variety of different companion radii. We find that our best-fit transit light curve has A% 2 ~-5 
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relative to a constant flux fit. Based on analysis of the noise properties of the KELT light curve and 
the number of trials we performed searching for a transit, we estimate that A^ 2 <-16 is generally 
indicative of a reliable detection, and thus this improvement is not significant. 

We then determine the fraction of trials that lead to a A% 2 greater than some threshold value. 
The results for A% 2 =9, 16, and 25 are shown in Fig. [51 We find that ~95% of MCMC realizations 
of transit models for companion radii >1.2R Jup lead to fits to our light curve that are excluded by 
our data, in the sense of producing a A% 2 that is worse by more than 16 relative to a constant fit. 
Therefore, we exclude with ~95% confidence that the companion transits if it has a radius larger 
than ~ l.2Rj up , and with ~75% confidence if it has a radius larger ~Rj up . We conclude that while 
transits of a Jupiter-radius companion are unlikely, they are not definitively excluded. 



6. STELLAR PARAMETERS 

We have made multiple determinations of the stellar parameters of the host star, using sev- 
eral different sets of data and analysis methods, described below. The results are summarized in 
Table [6l We note that, although the different determinations are generally mutually consistent, the 
uncertainties associated with each are simply formal statistical uncertainties, which have not been 
externally calibrated. We expect that these formal uncertainties are likely underestimates of the 
true uncertainties. Therefore, we conservatively choose to report the median of the three highest 
resolution spectroscopic results as our best estimate of the stellar parameters, and take the uncer- 
tainty as the standard deviation of the three estimates. The final stellar parameters we adopt are 
effective temperature r e ff = 6186 ± 92 K, surface gravity logg = 3.89 ± 0.07 (cgs), and metallicity 
[Fe/H] = -0. 15 ± 0.04. These and other properties of the star are listed in Tabled 



6.1. Fitting of Spectral Lines 



We analyzed the extracted APO 3.5-m spectra to det ermine the stellar p roperties in a careful 



hand-guided analysis according to the technique s used by 



Laws et al 



(2003), which are described 



more fully (excepting recent improvem ents) in IGonzalez & Vanturd (11998b. Briefly, we make 
use of the line analysis code MOOG (|Snedenl Il973l updated version), the iKuruczl (| 19931) LTE 
plane-parallel model atmospheres, and equivalent width (EW) measurements of 62 Fe I and 10 Fe 
II lines to determine the atmospheric parameters r eff , lo gg, microturbulence and [Fe/H]. The 
formal uncertainties were calculated using the method in IGonzalez & Vanturd (I1998|) . The values 
are listed in Tabled 



As a check, we performed a second analysis of the APO spectra using the code Spectroscopy 
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Made Easy (SME; see IValenti & Fischer! |2005|) . SME is an IDL-based program that uses syn- 
thetic spectra and least-squares minimization to determine the stellar parameters (e.g., r eff , log g, 
[Fe/H], vsinz, etc.) that best fit an observed spectrum. To constrain the stellar parameters, we 
analyzed three wavel ength regions (5 160 -5 190 A, 6000 -6200 A, and 6540 -6590 A) used by 
Stempels et al. d2007h . The first region is sensitive to log g. The second region contains a large 
number of spectral features of different elements, and is sensitive to [M/H] and v sin i. The third re- 
gion contains Ha, and the broadening of the outer wings of this line is sensitive to T eS . We fitted all 
three regions simultaneously using SME to estimate the stellar parameters of TYC 1240-00945-1. 
SME was unable to determine v sin i to a level finer than the velocity resolution of the APO 3.5-m 
spectra (~ 9 km s" 1 at R ~ 31500). We derive d parameters that agreed with those determined from 
the same spectra using the Laws et al.1 (|2003|) methodology. The values are listed in Table [61 



The stellar parameters were also verified using the ESO 2.2-m FEROS spe ctra. Measurement s 
of the equivalent widths were carried out automatically using the ARES code (ISousa et al.ll2007|) . 
Given the high S/N and broad spectral range of the spectrum, results were obtained for a large 
number of atomic li nes. However, after a careful inspection, only 21 Fe I and 9 Fe II lines (from 



the list in Table 2 of Ghezzi et al. 



201(h ) were considered sufficiently relia ble to be used in t he de 
termination of the stellar parameters. Applying the technique described in iGhezzi et al.l (|2010h . the 
following results were obtained: T eS = 6186 ± 82 K, logg = 4.01 ± 0. 17, f, = 1. 26±0.17kms -1 , and 
[Fe/Hl = -0.14 ±0.08, where the formal uncertainties were calculated as in Gonzalez & Vanture 
dl998b . 



The projected rotational velocity of TYC 1240-00945-1 was esti mated from the hig h-resolution 
FEROS spectrum using a technique similar to the one described in iGhezzi et al.l (120091) . The ex- 
pectation from FEROS simulations is that the high oversampling of the line spread function for 
the FEROS spectrum allows us to probe to much lower v sin i than achievable with the APO 3.5-m 
spectra, even though the FEROS resolving power is only moderately higher. We measure v sin i 
by simultaneously fitting the macro-turbulence velocity and v sin i for three moderately strong Fe 
I spectral lines. A grid of synthetic spectra was generated, varying v sin i, the macro-turbulence 
velocities and the adopted [Fe/H], the latter by 0.05 dex around the mean value given above. Small 
adjustments in the continuum level under 0.4% were allowed, to account for possible errors in 
the normalization process. In addition, small shifts in the central wavelengths of the Fe I lines 
were needed in order to properly match the observed lines. Values for v sin i and macro-turbulence 
were determined separately for each of the Fe I lines considered, based on standard reduced-% 2 
minimization. The results obtained for the three Fe I lines were consistent, yielding a v sin i in the 
range 1.1 -3.2km s~\ and macro-turbulence in the range 4. 5-4. 7km s" 1 . The l atter values are in 
good a greement with the macroturbulence velocity derived from Equation 1 in IValenti & Fischer 
(|2005|) and T e s = 6186 K. Our best estimate for v sinz was computed as the mean of the three 
values, yielding v sin i = 2.2 ± 1.5 km s" 1 , where the uncertainty is the RMS value; this RMS scat- 
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ter is approximately equal to the intrinsic uncertainty of the fitting procedure, which is typically 
1-2 km s" 1 . Note that when we tried recovering v sin i from simulations of FEROS spectra at 
the S/N of the TYC 1240-00945-1 spectrum, we found that even lower v sinz would indeed be 
detectable at the FEROS resolution. However, as discussed below in §17.2[ the lower la limit of 
0.7 km s" 1 leads to a very long rotation period which is astrophysically unlikely; it is more probable 
that the true v sin i lies within the upper half of the estimated range from the fit. 

We searched the FEROS spectra for any indication of spectral features from a secondary star 
blended with the primary, as might be expected if the RV signal were in fact caused by a nearly 
pole-on orbit of a low-mass stellar companion. To make a quantitative search for extra flux, we 
computed the difference of the normalized spectrum of TYC 1240-00945-1 with a tem plate FEROS 



spect r um of the primary star of the binary H P 20010, a well-studied F subgiant (|Balachandran 



19901 : iSantos et a l. 2004; iLuck & Heiterl 120051) with stellar parameters similar to those we derived 
for TYC 1240-00945-1. We examined the 8570 -8630 A region, where the spectrum has good 
continuum level determination, several spectral features (mostly due to Fe I), and where the con- 
trast ratio between an M dwarf and the primary would be relatively high, before the red-end fall-off 
in detection efficiency of the FEROS spectra (in this range, the S/N p er pixe l of T YC 1240-00945- 



1 and HD 20010 were high: 180 and 390, respectively). From the IPicklesI (| 19981) low resolution 



spectral library, we computed the expected ratio of fluxes between F8IV and MOV stars over this 
wavelength range to be 1.5%. We would expect that ratio to manifest as a difference in line ratios 
between the template and target spectra, with the M dwarf's flux filling up the cores of the F star's 
lines. However, the difference spectrum shows no detectable systematic offsets at the locations of 
HD 20010's lines; rather, the difference is evenly distributed around zero, with a standard devia- 
tion of 1.0%. The difference spectrum is shown in Fig. [6l This amount of deviation is expected 
since there is uncertainty in picking a template which would exactly match TYC 1240-00945-1. 
Thus there is no evidence for an MOV contaminating spectrum, although much cooler M dwarfs 
would provide less than 1.5% contaminating flux and might not be visible given the noise in our 
measurement. 



6.2. Spectral Energy Distribution Fitting 



As an additional check on the parameters of TYC 1240-00945-1, we performed a model 
atmosphere fit to the observed spectral energy distribution (SEP) f rom the optical fluxes from 
HAO ( ffT2l) and near-IR fluxes from 2MASS Jskrutskie et alJl2006l) . The absolute photometric 
measurements in the g'r'i'JHK s passbands (see Table© were converted to physical fluxes using the 
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published SDSS0 and 2MASS0 zero- points, together with p ublished color-dependent corrections 
to the passband effective wavelengths ( Moro & MunarilboOO ). The model atmospheres used in the 
fitting are the NextGen atmospheres of lHauschildt et al.1 (|1999|) . which are gridded in r eff by 100 K, 
in logg by 0.5 dex, and in [Fe/H] by 0.5 dex. We performed a least-squares fit of this model grid to 
the six flux measurements, with the extinction Ay and the overall flux normalization as additional 
free parameters. 

We initially allowed all of the variables- T eS , logg, [Fe/H], A v , and flux normalization- to be 
fit as free parameters. We limited the A v to a maximum of 0.65, corresponding to the maximum 
line-of-sight extinction as determined from the dust maps of ISchlegel et al.l (|1998[) . The resulting 



fit is shown in Fig. with T eS = 6400!^° K, A v = 0.6^ f 5 , logg = 3.5 ± 1 .5, and [Fe/H] = 0.0 ± 2.0. 

These values are consistent with those derived spectroscopically. However, the available pho- 
tometry does not strongly constrain the stellar parameters as there is a very strong degeneracy in 
the SED fit between T s q and Ay, due to the lack of absolute flux measurements at wavelengths 
bluer than 0.5 ftm. Thus, we re-fit the fluxes with r e ff fixed at the spectroscopic value of 6186 K, 
[Fe/H] fixed at 0.0, and logg fixed at 4.0; the only remaining free parameters are A v and the nor- 
malization. In this way we use the photometry to strongly constrain the line-of-sight extinction. 
The resulting best fit, with A v = 0.40 ± 0.05, is displayed in Fig. [71 



Adopting this A v , which implies E(B-V) = 0.13 using the reddening law of iBessell & Brett 



( 198 8]), we can check r e ff from the broadband colors alone, using the recent color calibrations of 
Casagrande et al.l (|2010|) . For example, from the J-K s color we find r eff = 6147 K, while from the 
V -K s color we obtain T e ^ = 6299 K. Thus, given a reasonable estimate of A v , even when we use 
individual colors instead of fitting them all simultaneously, the T e & estimates are consistent with 
the spectroscopically determined value to within ~ 100 K. 



6.3. Final Determination of the Stellar Parameters and Companion Parameters 



We determine the mass and radius of the parent s tar, TYC 124 0-00945-1, from T eff , logg, 
and [Fe/H] using the empirical polynomial relations of iTorres et al. (I2010I) . which were derived 
from a sample of eclipsing binaries with precisely measured masses and radii. We estimate the 
uncertainties in M* and R* by propag ating the uncertaint ies in T eS , logg, and [Fe/H] (see Table [2]) 
using the covariance matrices of the Torres et al. ( 201p|) r elations kindly provided by G. Torres. 



Also, since the polynomial relations of 



Torres et al 



(|2010|) were derived empirically, the relations 



^http : //www . sdss . org/dr 7 /algorithms /f luxcal . html| 

^http ://ssc.spitzer.caltech. edu/ documents/ cookbook /html /cookbook- node 2 07 . html| 
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were subject to some intrinsic scatter, which we add in quadrature to the uncertainties propagated 
from the stellar parameter measurements. The final stellar mass and radius values we obtain in this 
way are M* = 1 .37 ± 0. 1 1 M Q and R* = 2.2O!°;|fl . 

Using the derived value of M*, we estimate a minimum mass (i.e., for sin i = 1 where i is 
the orbital inclination) for the companion, MARVELS-lb, of m m j n = 28.0 ± \ .5M Jup , where the 
uncertainty is dominated by the uncertainty in the primary mass. In fact, the mass function, 

(msim) 3 9 , /9 , 

is more precisely determined. We find (msinz) 3 /(M* + m) 2 = (9.75 ±0.32) x 1O~ 6 M . With our 
adopted value of M*, we can also estimate the semimajor axis a = 0.071 ± 0.002 AU, assuming an 
edge-on orbit; for less inclined orbits, the semimajor axis is larger. 

The small minimum mass of the companion positions it as a good short-period brown dwarf 
desert candidate. In order for it to be a low-mass star rather than a brown dwarf, the orbital 
inclination would have to be close to face-on. In order to explore further the probability that the 
companion has a mass greater than the hydrogen burning limit, we conducted a Bayesian analysis 
to estimate the posterior probability dist r ibutio n for the companion mass, using the methodology 



described in Section 7 of Fleming et al.l (120 1 QQ : an MCMC chain is constructed starting from a 
distribution of stellar parameters an d error bars as ado pted for TYC 1240-00945-1 in Table [2l 
stellar masses are determined using ITorres et al. (boioh . and companion masses are determined 



using a random distribution of inclinations. This analysis assumes a uniform distribution in cos i, 
includes uncertainties on the orbital and host star parameters, and adopts priors on the luminosity 
ratio and mass ratio for the companion. 

Of course, the posterior distribution of the true companion mass d epends on our adopted 



prior for the companion mass ratio distribution (e.g., Iho & Turnerlboioh . Given that few brown 
dwarf companions are known, the constraints on the companion mass ratio distribution in the mass 
regime of interest are poor. Indeed, this is what makes this object interesting, and this distribution 
is precisely what we would like to infer from a larger ensemble of similar detections. Neverthe- 
less, we can adopt various simple and plausible forms for the mass ratio distribution, and then 
use these to infer posterior probability distributions for the true mass. From Doppler surveys for 
exoplanets, it is known that Jupiter-mass companions are significantly more common than brown 
dwarf companions, and that the frequency of planetary companions declines for larger masses, 
such that the mass fu nction is roughly uniform in the logarithm of the planet mass for m < 10 Mj 



(ICumming et al.ll2008|) . It is not known if this form holds for companions with mass significantly 
larger than ~ 10 Mj, but it is clear that the frequency of companions in the brown dwarf regime 
must reach a minimum at some point and then rise again, given that M dwarf companions with 
masses just above the hydrogen burning limit are known to be more common than brown dwarf 
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companions. IGrether & Lineweaverl (|2006[) found that this minimum (the driest part of the brown 
dwarf desert) occurs at a companion mass of 31+jf Mj. Thus the minimum mass of MARVELS- 
lb is near the minimum of the companion mass function, and prior mass ratio distributions that 
are falling, flat, or perhaps rising shallowly in logg are all equally plausible (see Figure 11 of 
Grether & Lineweaver 2006). 



We therefore consider five different priors on the companion mass ratio distribution: dN / d log q oc 
q" 1 , oc logg, constant, oc q, and oc q 2 . The first thre e are falling or constan t with l ogg, and the lat- 
ter two are rising with logg. From the results of IGrether & Lineweaver! (|2006), we believe the 
first three are the most plausible, while the first four almost certainly bracket the likely range of 
distributions for companions close to the relevant regime. The resulting cumulative probabilities 
for the companion mass for the five different priors are plotted in Figure [8] For the three favored 
priors, we conclude that at > 90% confidence the actual mass is below the hydrogen-burning limit. 
For the prior that is uniform in (linear) mass ratio, d/V/dlogg oc q, there is a ~ 25% probability 
that the companion is in fact a low-mass star, whereas it is only for the assumption of relatively 
steeply rising mass ratio distribution (dN/dlogq oc q 2 ) that the companion is more likely to be a 
star. Again, we do not believe such a distribution is very likely to be correct for this regime of 
companion mass, but given the poor constraints, we cannot absolutely exclude it either. Finally, 
we note that for the last two priors, the precise form of the posterior distribution depends on our 
imposed constraint on the luminosity ratio, which is somewhat uncertain. 



With a reddening of E(B — V) = 0.13 ( §16.21 ), the system is evidently seen much of the way 
through the full reddening along this line of sight, which from the lSchlegel et al.l (|1998|) dust maps 
is E(B-V) = 0.186. The physical distance of the system can be estimated from its luminosity 
and apparent magnitude. First we compute the bolometric magnitude of the star as M bo i = 4.74- 
2.51og(L/L Q ), where 4.74 is the bolometric magnitude of the Sun. The luminosity is calculated 
from the Stefan-Boltzmann law applied to the r eff and st ellar radius calculated above , and we adopt 
a BC V =-0.17 as appropriate for its spectral type (e.g. jKenyon & Hartmannlll995|) . The absolute 
magnitude is therefore 2.91. Adopting Ay = 0.4 ±0.05 ( §16.21) . this yields a distance d = 280 ±30 pc. 



6.4. Expected Stellar RV Jitter 



Starspots and motions of the stellar surface are possible astrophysical sources of noise that 
can interfere with searches for compa nion RV signals. These so u rces a r e commonly ref e rred t o 
as " jitter", and are exp l ored b y, e.g., Isaar et al.1 d 1998b . Iwrightl hp05\) . Lagrange et aD d2009h . 
and llsaacson & Fischer] (|2010|) . For late F dwarfs of B- V > 0.5, they find typical jitters in the 
~10m s" 1 range, with the most extreme outliers at ~100m s" 1 . 
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TYC 1240-00945-1 is slightly evolved, so one wonders whether it might experience larger 
jitter than for F dwarfs. However, it still lies at B—V and M v below and redward of th e instability 
strip (for a review of the position of the strip, see, e.g.. ISandage & Tammannl 120061) . and shows 
no signs of activity based on the time-series photometry (Section [5])- so one shouldn't expect 
multi-periodic pulsations at the level o f, e.g., the ~400m s " 1 RV jitter of the brown dwarf-hosting, 
instability strip, A9V star HD 180777 (|Galland et al.ll2006l) . Rather, F stars with stellar parameters 
similar to TYC 1240-00945-1 can be fairly quiet in ter ms of the RMS sc atter attributable to RV 
jitter: ~4-5 m s _1 in the case of the F6 star HP 605 32 dDesort et alJbood) . and ~10m s" 1 in the 
case of the F7 star HD 89744 ( Korzennik etaPboOO ). 



We conclude that the levels of RV jitter expected for this combination of stellar parameters are 
too low to be responsible for the K = 2.533 ±0.025 km s" 1 of the TYC 1240-00945-1 RV signal, 
although they could be a contributor to the extra error we have regarded as systematics in the RV 
analysis. 



7. DISCUSSION 



7.1. Evolutionary state of the host star 



In Fig. [9] we compare the spectroscopically measured T e ff and logg of TYC 1240-00945-1 
(red erro r bars) against a theore tical stellar evolutionary track from the Yonsei-Yale ("Y 2 ") model 
grid (see IDemarque et al.ll2004l and references therein). The solid curve represents the evolution 
of a single st ar of mass 1.37M m (the mass of TYC 1240-00945-1 inferred from the empirical 
calibration of iTorres et alj|2010c see above) and metallicity of [Fe/H]=-0.15 (as determined spec- 
troscopically), starting from the zero-age main sequence (lower left corner), across the Hertzsprung 
gap, and to the base of the red-giant branch. Symbols indicate various time points along the track, 
with ages in Gyr labeled. The dashed curves represent the same evol utionary track but fo r masses 
±0.1 1M Q , representative of the la uncertainty in the mass from the Torres et al.1 (|2010h relation. 
The filled gray region between the mass tracks therefore represents the expected location of a star 
of TYC 1240-00945- l's mass and metallicity as it evolves off the main sequence. We emphasize 
that we have not directly measured the mass of TYC 1240-00945-1, and thus we are not attempt- 
ing to test the accuracy of the stellar evolutionary tracks. Rather, our goal is to use these tracks to 
constrain the evolutionary status of the TYC 1240-00945-1 system. 

The spectroscopically measured T eS , logg, and [Fe/H] place TYC 1240-00945-1 near the 
beginning of the subgiant phase, just prior to crossing the Hertzsprung gap to the base of the red 
giant branch, with an estimated age of ~ 3 Gyr. 



We can also take advantage of the information provided by the MARVELS input catalog to 
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place the host star on a n RPM diagram, taking c olors from the 2MASS catalog, and proper motions 
from the GSC2.3 (see iGould & Morganll2003l for an example of how RPM can be used to help 
differentiate giants from dwarfs). In Fig. \\0[ we show that the /-band RPM (RPMj = J + 5 log fi) 
is most consistent with the host star being a dwarf or subgiant, as it falls well away from the region 
of the RPM diagram dominated by giant stars. 



7.2. Tidal Effects 



Given the relatively large mass ratio and short period of the TYC 1240-00945-1 system, tidal 
interactions between the star and MARVELS- lb could be important- given the roughly ~3 Gyr age 
of the host star, is the system likely to be tidally synchroni zed? We follow exactly the same analysis 
of the tidal interaction as detailed in Fleming et al.1 (|2010|) . which uses the tidal quality factor of 
the star, Q'^, as a free parameter in the equations for the decay of the companion's semimajor axis 
over time and the relation of t he primary's rotatio nal frequency to the companion's orbital angular 
momentum (Eqs. 5 and 6 in Fleming et al.1 |2010|) : together, the equations permit a solution for 
the amount of time required for tidal sync hronization. Note t hat if the primary's rotation never 



sync hronizes, the two bodie s may merge (ICounselmanl 1 19731 : iLevrard et al.l l2009i : Uackson et al 



2009J). As in Fleming et al.l (I2010|) . we have examined the tidal evolution of this system in the 
range 10 4 < <2* < 10 10 , for a range of values of the inclination of the secondary's orbit to the line 
of sight from i = 0° (face-on) to z = 90° (edge-on), adjusting the mass and rotation period using the 
measured values of v sin i and R* from §(6] We set the primary's equator to be in the same plane as 
the secondary's orbit, but this decision does not affect our results. 



In Fig.QT]we show the synchronization and merging times from Eqs. 5 and 6 of lFleming et al. 



(120101) . over the <2* and i parameter space defined above. The curves are isochrones in the (<2*, i) 
parameter space, so if the TYC 1240-00945-1 system has a (<2*, i) combination that lies above a 
given isochrone T sync / merge , then the system will take longer than T sync / merge to synchronize or merge. 
Isochrones are plotted for T sync / merge = 0.01, 0.1, 1, and 10 Gyr. 

We consider three models: the best-fit stellar parameters (solid curves); one in which v sin i = 
3.7km s" 1 , M* = 1.43M©, and i?* = 2.44 R Q (dotted curves); and one with v sin i = 0.7km s" 1 , 
M* = 1.32M , and R* = 2.00 R Q (dashed curves). The latter two cases represent models where the 
parameter sets were adjusted in opposite directions in an attempt to have the two models span a 
maximal amount of (£)*, i) parameter space, while still maintaining the parameters within the un- 
certainties. Thus, the uncertainty on the four synchronization/merging isochrones is approximately 
indicated by the region between the dotted and dashed lines (though it is not a perfect indication 
of the multi-parameter uncertainty envelope, as is evident from the fact that the dotted and dashed 
lines cross). 
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Note if one makes a trial assumption for the value of the inclination i, then given our mea- 
surement of v sin i, one may infer the true rotational velocity v of the stellar surface. At some small 
inclination, close to a face-on orbit, this will yield a v so high that the primary's rotational fre- 
quency is already spun up to tidal synchronization (and for the improbable case of an inclination 
even smaller than this, the primary's rotational frequency is higher than the secondary's orbital 
frequency, a scenario we do not explore here, but which would result in gradual spindown of the 
primary's rotational frequency until it matched with the orbital frequency of the secondary). For 
each case of v sin i that we investigated, the value of the inclination which corresponds to present- 
day tidal synchronization is visible on Fig.[TT|as a vertical asymptote towards which the isochrones 
converge. For inclinations closer to edge-on, the secondary still is in the process of spinning up the 
primary. 

Next consider the best fit (solid curves) and maximum v sin i (dotted curves) cases. We find 
that for a wide range of low (<2*, i) combinations, the secondary quickly spins the primary up 
to synchronization in less time than the ~ 3 Gyr age of the host star. However, this alone, while 
suggestive, is not conclusive proof that such a synchronization has occurred. As this is an evolved 
F star, the radius has recently expanded, complicating any interpretations of the system's history. 
Furthermore, for <2* ~ 10 7 , the synchronization time is about the age of the system. 

For the minimum v sin i cases (dashed curves), the rotational period of the star is very large, 
~ 150 days. While this period may not be physical, it is formally permitted by the observations. 
With such slow rotation, the companion may merge with the star before synchronization is finished. 
This possibility of the synchronization timescale exceeding the merging timescale occurs when 
i > 54° (note there is no feature in Fig. [TT]at the z = 54° transition, because in our simplified model 
a companion can reach the stellar surface and synchronize the star's rotation period, or move just 
inside the surface and merge) . Undoubtedly the behavior of such a compact system is not well- 



modeled by Eqs. 5 and 6 of [Fleming et al.l (|2010|) , but we cannot rule out the possibility that 



MARVELS- lb will eventually merge with the host star. 



8. SUMMARY 

In a search through the first year of SDSS-III MARVELS data, we have discovered MARVELS- 
lb, a candidate brown dwarf companion to the V ~ 10.6 star TYC 1240-00945-1 with a velocity 
semiamplitude of K = 2.533 ±0.025 km s" 1 and an unusually short period of 5.8953 ± 0.0004 d. 
Radial velocity data from several observatories confirm the Doppler variability, and high-resolution 
spectroscopic observations indicate that the host is a mildly evolved, slightly subsolar metallicity 
F star with T eff = 6186 ±92 K, logg = 3.89 ±0.07, and [Fe/H]=-0. 15 ±0.04, with an inferred mass 
of M* = 1.37 ±0.11 M . The minimum mass of MARVELS- lb is 28.0 ± l.5M Jup , implying that 
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it is most likely in the brown dwarf regime. We see no evidence for spectral lines from the com- 
panion in the high-resolution spectra, implying that the companion is not an M dwarf with an orbit 
extremely close to pole-on. Comprehensive, precise relative photometry indicates no variability at 
a level of > 1% on time scales of hours to years. Phasing to the period of MARVELS- lb as well 
as the first harmonic, we can place an upper limit on the amplitude of coherent photometric vari- 
ability of ~0.05%. Under many (but not all) of the potential combinations of system parameters, 
this short-period system is likely to have tidally synchronized, given the estimated ~ 3 Gyr age of 
the host star. 

The a priori transit probability of MARVELS- lb is quite high, ~ 14%. Although we find no 
evidence for transits, we also cannot definitively rule them out for likely MARVELS- lb radii of 
r^Rjup- The transit ephemeris is T c = 2454936.555 ± 0.024 (BJD TDB ), with an expected transit 
depth of r^0.2%(r/R Jup ) 2 , and a duration of ~6.5 hours for a central transit. 

We believe this candidate highlights the great promise of MARVELS as a factory for finding 
the rare companions that populate the brown dwarf desert. The primary goal of the MARVELS 
survey is to monitor ~ 10 4 main sequence and subgiant stars with velocity precision sufficient to 
detect Jovian companions with periods of less than a few years. As such, MARVELS is uniquely 
and exquisitely sensitive to massive but rare companions. MARVELS- lb is the first of a number of 
brown dwarf candidates we have identified in the MARVELS data obtained to date, and we expect 
to uncover several additional such systems as the survey progresses. 

Funding for the MARVELS multi-object Doppler instrument was provided by the W.M. Keck 
Foundation and NSF with grant AST-0705139. The MARVELS survey was partially funded by 
the SDSS-III consortium, NSF grant AST-0705139, NASA with grant NNX07AP14G and the 
University of Florida. Funding for SDSS-III has been provided by the Alfred P. Sloan Founda- 
tion, the Participating Institutions, the National Science Foundation, and the U.S. Department of 
Energy. The SDSS-III web site is http : / /www . sdss3 . org/| SDSS-III is managed by the As- 
trophysical Research Consortium for the Participating Institutions of the SDSS-III Collaboration 
including the University of Arizona, the Brazilian Participation Group, University of Cambridge, 
University of Florida, the French Participation Group, the German Participation Group, the Michi- 
gan State/Notre Dame/JINA Participation Group, Johns Hopkins University, Lawrence Berkeley 
National Laboratory, Max Planck Institute for Astrophysics, New Mexico State University, New 
York University, the Ohio State University, University of Portsmouth, Princeton University, Uni- 
versity of Tokyo, the University of Utah, Vanderbilt University, University of Virginia, University 
of Washington and Yale University. The Hobby-Eberly Telescope (HET) is a joint project of the 
University of Texas at Austin, the Pennsylvania State University, Stanford University, Ludwig- 
Maximilians-Universitat Miinchen, and Georg-August-Universitat Gottingen. The HET is named 
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Table 1: SDSS-III MARVELS Radial Velocities for TYC 1240-00945-1 
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Table 2: TYC 1240-00945-1: Parameters of the Star 

Parameter Value 



Spectral Type 

g' 
r' 

i' 
B 
V 
R c 

lc 

J 2MASS 
H-2MASS 
KjMASS 
Teff 

[Fe/H] 

Mass 

Radius 

Ay 

Distance 
v sin i 



F9IV-V 
10.821 ±0.013 
10.436 ± 0.007 
10.324 ± 0.013 
11.230 ±0.025 
10.612 ±0.025 
10.242 ±0.011 a 
9.916 ±0.011 a 
9.395 ± 0.018 
9.112 ±0.016 
9.032 ± 0.017 

6186 ± 92 K 
3.89 ± 0.07 (cgs) 

-0.15 ±0.04 
1.37 ±0.11 M Q 

2-2O^[];22^0 
0.40 ±0.05 
280 ±30 pc 
2.2 ±1.5 km s" 1 



a R c I c are transformed magnitudes based on g'r'i', using the transformation equations of 



Smith et al. 



(2002). 
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Table3: HET Radial Velocities for TYC 1240-00945-1 
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Table 4: SMARTS Absolute Radial Velocities for TYC 1240-00945-1 



BJD TDB 


Absolute RV 
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Table 5: MARVELS- lb: Parameters of the Companion 



Parameter 



Value 



Minimum Mass 
a 
K 
P 

^prediction for transit 

ecosw 
esinu; 



28.0 ± 1.5 M Jup 

0.071 ±0.002 AU 

2.533 ± 0.025 km s" 1 

5.8953 ± 0.0004 d 

2454936.555 ± 0.024 (BJD TDB ) 

u.uu _ 010 

-0 003 +() 0()8 
w.wwj _ ooog 



Table 6: TYC 1240-00945-1 Individual Determinations of Stellar Parameters 
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SED fit to photometry 
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TYC 1240-00945-1 MARVELS orbitol fit 
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TYC 1240-00945-1: MARVELS residuols after subtracting orbital fit 

u > 



i 



I 



HIP 14810: MARVELS residuals for nearby RV reference star 
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Fig. 1.— Top: MARVELS RV data and Keplerian orbital solution for TYC 1240-00945-1. Beam 
1 is shown with blue filled squares, and beam 2 with green open squares. Center: The residuals 
for TYC 1240-00945-1, equal to the RVs from the top panel minus the orbital fit. Bottom: The 
residuals for HIP 14810, a star with a known two-planet RV signal, observed through a nearby 
fiber during the same exposures as those plotted for TYC 1240-00945-1. 
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Fig. 2. — Phase-folded Keplerian orbital solution and RV residuals for TYC 1240-00945-1. Blue 
squares and green squares are MARVELS discovery data, red circles are HET data, and purple 
triangles are SMARTS data. Error bars have been scaled up by the methodology in §14.11 The 
bottom panel shows the residuals between the data points and the orbital solution. Note that the 
HET and SMARTS data were not used in the Keplerian fit, and so provide an independent check 
of the quality of the MARVELS data. 
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Fig. 3.— Top: KELT North light curve for TYC 1240-00945-1. Bottom: Lomb-Scargle peri- 
odogram of the KELT data, showing no evidence for any significant periodicities for periods of 
P= 1 — 10 days, including the period of MARVELS- lb (vertical dashed line) and the first harmonic 
(vertical dotted line). 



-33- 



-0.02 



0) 

cd 
S 

> 
•i— i 

+-> 








0.2 



0.4 0.6 
Phase 



0.8 



Fig. 4.— Grey points: the KELT light curve for TYC 1240-00945-1, phased to the period of 
MARVELS-lb (5.8953 days). Black points: the phased KELT light curve, binned using bin size 
A0 = O.O4. 
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Fig. 5. — Probability that transits of MARVELS- lb are excluded at levels of A% 2 = 9 (solid red), 
16 (dotted blue), and 25 (dashed green), based on the analysis of the KELT photometric dataset, 
as a function of the radius of MARVELS- lb. Also shown is the case for A\ 2 = 16, but assuming a 
box-shaped transit (black, long dashes) instead of a limb-darkened light curve model. 
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Fig. 6. — The normalized spectrum of TYC 1240-00945-1 minus the normalized spectrum of 
HD 20010 (an F8IV star with similar stellar parameters), in the wavelength range 8592-8626A. 
Locations of some spectral lines of a late F subgiant are indicated by red dots. 
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Fig. 7. — Top: NextGen model atmospheres SED fit to the color photometry of TYC 1240-00945- 
1, allowing all parameters to vary. Bottom: SED fit to the color photometry of TYC 1240-00945-1, 
with only Ay and the normalization as variables (r eff , [Fe/H], and logg locked). 
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Fig. 8. — Cumulative probability that the mass of MARVELS-lb is less than a given mass, in 
units of solar masses. These probabilities account for the uncertainties and co variances between 
the parameters of the Keplerian orbital fit, the uncertainty in the host star mass, the assumption 
of a uniform distribution of cos i, and the adoption of five different priors for the distribution of 
companion mass ratios dN/d\ogq. 
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Fig. 9. — The evolutionary track for an object with M = 1.37±O.11M , at [Fe/H]=-0.15. Ages of 
1.0, 2.0, 3.0, 3.2, 3.3, and 3.5 Gyr are indicated as dots. The possible tracks for up to a la deviation 
in the mass are shown by the shaded region. The stellar parameters for TYC 1240-00945-1, with 
lcr error bars, are shown by the cross. 
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Fig. 10. — ./-band reduced proper motion versus J — H color. Stars from the RAVE DR2 
(jZwitter et al.ll2008r) with galactic latitude, 20° < \b \ < 30°, and with measured spectroscopic prop- 
erties are shown. The RAVE stars are color coded by luminosity class such that giants (logg<3.5) 
are red, dwarfs (logg > 4. 1) are green, and subgiants (4 . 1 >log g > 3.5) are blue. The polynomial 
relation (solid line) defined from [Collier Cameron et al.l (I2007|) discriminates the dwarf star popu- 
lation from the giant star population in this plane. TYC 1240-00945-1, plotted as the large black 
circle, is consistent with being a dwarf or subgiant. 



90 60 30 
Inclination (°) 

Fig. 11. — Contours of the time (in Gyr) to synchronize the primary's rotational period to the 
orbital period, or for the companion to merge with its host star. Solid curves correspond to the best 
fit, dotted curves to the case with v sin i, M*, and R* each set at the tops of their la uncertainty 
ranges, and dashed curves to the case with v sin i, A/*, and R* each set at the bottoms of their 
la uncertainty ranges. Merging is only possible in the latter case when i> 54°; therefore those 
portions of the dashed curves represent the time to merge. 
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